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Abstract 

We describe systematic ranging, an orbit determination techniqne especially 
snitable to assess the near-term Earth impact hazard posed by newly dis¬ 
covered asteroids. For these late warning cases, the time interval covered by 
the observations is generally short, perhaps a few honrs or even less, which 
leads to severe degeneracies in the orbit estimation process. The systematic 
ranging approach gets aronnd these degeneracies by performing a raster scan 
in the poorly-constrained space of topocentric range and range rate, while 
the plane of sky position and motion are directly tied to the recorded ob¬ 
servations. This scan allows ns to identify regions corresponding to collision 
solntions, as well as potential impact times and locations. From the probabil¬ 
ity distribntion of the observation errors, we obtain a probability distribntion 
in the orbital space and then estimate the probability of an Earth impact. 
We show how this techniqne is effective for a nnmber of examples, inclnding 
2008 TC 3 and 2014 AA, the only two asteroids to date discovered prior to 
impact. 
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1. Introduction 


When an asteroid is hrst discovered, the few available observations rarely 
allow the determination of its orbit. Surveys usually capture 3-5 repeat 
images of the same area of sky with a typical return time of 15-30 minutes. 
Over this time interval, an object in the held of view moves with respect to 


background stars at a roug 


ily linear rat e . Th us, the short arc of available 


data, also called a tracklet flKubica et al.l. 1200711 . provides an estimate of the 


angular position and motion of the object i n the plane of sky, which i s not 
enough to determine a six-parameter orbit (iMilani and Knezevid 1200511 . 

Despite the degeneracies in the orbit estimation process, it is important 
to recognize potentially hazardous objects shortly after discovery. In fact, the 
only two objects discovered prior to an Earth i mpact, namely 2008 TC.q and 


2014 AA, were both discove red by R. Kowalski flKowalski et al 


2008 


201411 


of the Catalina Sky Survey (ILarson et ahl. Il998l) only about 20 hours before 
striking the Earth. On one hand, 2008 TC 3 was quickly recognized as poten¬ 
tial impactor and so it was extensively observed before the impact. The ac¬ 
quired observations allowed the estimation of the orbit of 2008 TC 3 as well as 


some physical c 


Scheirich et al. 


larac terization fjjenniskens et al. 


2009 


Kozubal et al 


2011 


201011 . Using the available data, it was poss ible to predict the 
impact that took place above Sudan on 2008 October f Jenniskens et al. 


^ http: //neo. jpl.nasa.gov / news /2008tc3.html 
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20091) . On the other hand, 2014 AA was not immediately recognized as a 
possible impactor and so only seven astrometric observations over about 70 
minutes wer e obtained before the object fell into the Atlantic Ocean on 2014 


January 2 g flCheslev et al. 


20151 


w: 


2015 


len the standard diff e rentia l correction procedure (e.g., 


Farnocchia et ah 


Milani and Gronchi l2010l) to find a least-squares orbit fails, other 


method s can be used to asses s the o rbital probability distribution. In par¬ 
ticular, iMuinonen and Bowelll (1199311 show how to put a probability density 
on the phase space of the orbital elements by using Bayesian inversion the¬ 
ory. By using a Monte Carlo approach one can sample the orbit phase space 
and thereby derive the probability distribution from the observation errors 
corresponding to the sampled orbit. 

The available observations directly constrain the position and motion of 
the asteroid in the sky while the distance between the asteroid and the ob¬ 
server (topocentric range) is poorly constrained. Thus, ranging methods are 
the prefer red Monte Carlo appr oach when only a short arc of observations is 


available. 


Virtanen et al. 


(1200 ll) describe a method called statistical ranging 


that allows one to generate orbital samples compatible with the observational 
data. After choosing two observations, the topocentric ranges at the epochs 
of the two observations are rando mly sampled and a corresponding orbit 


computed. 


Oszkiewicz et al. 


(1200911 rely on the statistical ranging technique 


and use Markov chains to generate an unbiased sequence of orbital samples 
distributed according to the probability distribution coming from Bayesian 


^http://www.nasa.gov/jpl/asteroid/first-2014-asteroid-20140102 
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inversion t 


Cheslev 


l eory. 


(120051) introduces a technique called systematic ranging, which 


in contrast to Monte Carlo techniques systematically explores a raster in 
the topocentric range and range-rate space. This technique provides a geo¬ 
metric description of the orbital elements as a function of range and range 
rate. Moreover, systematic ranging allows one to identify regions of the 
phase space filled with impact solutions and the corresponding impact times 
and locations. In this paper we present a detailed description of systematic 
ranging and show how to derive a probability distribution on the range and 
range-rate space, which is then mapped to the orbital element space where 
impact probability estimates can be derived. 


2. Systematic ranging 

Systematic ranging relies on the fact that a short arc of observations yields 
a direct estimate of the plane of sky position (right ascension a and declina¬ 
tion S) and motion (d and h). These four scalar quant ities can be assembled 
toget her in the so-called attributable A = (a, <5,0;, 5) flMilani and Knezeyid . 
20051) . The topocentric range p and p are only marginally constrained, if at 
all. If p and p were known, we would haye a full description of the aster¬ 
oid’s topocentric position and velocity in polar coordinates (a, 5, d, h, p, p), 
which can be easily converted to a Cartesian heliocentric state if the position 
and velocity of the observer are known. Note that, to account for light- 
time correction, we correct the epoch of the Cartesian state from that of the 


^This technique was actually introduced by Tholen and Whiteley in 2002, but the paper 
was never published. 


4 















attributable by a quantity p/c, where c is the speed of light. 

To explore the orbit phase space, systematic ranging scans a suitably 
dense grid in the (p, p) space. Such grid needs to be large enough to contain 


the reasonably possible orbital conhgur ations. In pa r 


contain the so-called admissible region flMilani et ah 


icular, the grid must 


20041) . i.e., the values 


of (p, p) leading to bounded heliocentric orbits. For each grid point we fix 
the values of p = p* and p = pj and hnd the best £t value of the attributable 


Aji that minimizes the cost function: 


Q = 


( 1 ) 


where u is the vector of Observed — Comput ed astrometric residuals and W 
is the weight matrix flFarnocchia et ahl. 120151) . 

The constrained best-fitting solution can easily be converted to an orbit, 
which is in turn propagated to find upcoming Earth encounters. Moreover, 
if the observations contain photometric measurements we also compute the 
absolute magnitude for each grid point. 

Figured] shows the application of systematic ranging to 2014 AA by using 
the seven astrometric observations obtained by the Catalina Sky Survey prior 
to impact. The orbital elements are shown as a function of (p, p). The dashed 
curve encloses bou nded orbits and corresponds to the admissible region of 


Milani et ah 


(120041) . The dash-dotted line is for grazing impacts: the region 


on the left of the curve contains impacting solutions. 


3. Probability distribution on the range &: range rate space 

The raster scan presented in the previous section gives a geometric de¬ 
scription of how the orbital conhguration depends on the value of topocentric 


5 











Topocentric range rate [au/d] Topocentric range rate [au/d] 


Eccentricity Perihelion distance [au] 






Figure 1: Eccentricity, perihelion distance, inclination, and absolute magnitude as a func¬ 
tion of range and range rate for asteroid 2014 AA. The dashed curve corresponds to 
parabolic orbits and the dash-dotted line to grazing impacts. 
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Table 1: Station-specific weighting rules for the most productive discovery and short-term 
follow-up stations. For all the other stations we use weights at 1". 

range and range rate. The next step is to assign a probability distribution 
to this space. 

As is common practice, given a set of optical observations we assume 
that the observation errors v are normally distributed according to a weight 
matrix W, i.e., their probability density f^rr is: 

ferrite) OC exp( —0.5 W I') . 

It is typical to use a diagonal weight matrix where the individual weights are 
1/cr^ (cr is the uncertainty in the astrometric positions). Table [1] shows the 
data weights we currently use for the most productive discovery and follow¬ 
up stations. The chosen weights are based on our experience and account for 
the fact that discovery observations can sometimes be problematic because 
they are not targeted on the object. 

In principle, it is possible to consider the photometric residuals together 
with the astrometric ones. However, photometric measurements have a much 
larger uncertainty than that of the astrometry and can be affected by signih- 
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cant biases fjjuric et all l2002[) . Moreover, in case of a magnitnde trend, it is 
hard to tell if the varying Inminosity is dne to a rapidly changing topocentric 
distance rather than the asteroid’s rotation. Therefore, we do not nse the 
information obtained from the photometri c residnals. _ 


According to Bayesian inversion theory fiMninonen and Boweli 


1993 1. the 


posterior probability density fnnction for (p, p) is: 

fpostiPy P) 0^ ferr {^{P. P)) fprior (P) P) ) 


where fprior is a prior distribntion on the (p, p) sp ace 


is to some extent arbitrary and far from trivial fjJavnes 


T he sele ction of /, 


prior 


19681) . and yet can 


signihcantly affect the posterior probability distribntion, especially for short 
observation arcs. 

Whatever the choice of fprior, we always add a crnde constraint from the 
popnlation model by setting fprior = 0 for hyperbolic orbits, flagging objects 
where the astrometric errors favor an nnbonnded orbit. 

3.1. Jeffreys ’ prior 


fprior (^P, P) det ^ ^ 


A mathematically sonnd choice is Jeffreys’ prior flGranvik et al.l. l2009l) : 

d(pj) 5 (^ 7 )) ■ 

Note that the partials are the total derivatives of ly with respect to (p, p), 
which means that they acconnt for the fact that the best-£t attribntable A 
changes as a fnnction of (p, p). 


''^It is often convenient to use a logarithmic scale for p to achieve a better resolution at 
small topocentric distances. In that case the probability density has to be multiplied by 

























Among other things, Jeffreys’ prior secures the invariance of the probabil¬ 
ity distribution when changing variables. However, since there is more sen¬ 
sitivity of the residuals for small topocentric distances, Jeffreys’ prior tends 
to favor orbital conhgurations where the object is close to the observer, e.g., 
see the examples of 2015 CV and 2015 BU 92 in Sec. Ei 


3.2. Prior based on population model 
As already discussed by IVirtanen et al. 


( 1200 ll) . the prior could be based 


on the known asteroid population. Even better, one could use an orbit and 
size population model where the known population is extrapolated to remove 
survey selection effects and to reach completeness up to a given apparent 
magnit ude. One such m odel is the Pan-STARRS Synthetic Solar System 


Model Gray et al. 


future flGranvik et al. 


20111 ). a nd new po pulation models a re expected in the 


2014 1. From the 


Gray et al. 


(1201 ll) model we derived 


a probability density fpop{q,e,i, H) for perihelion distance q, eccentricity e, 
inclination i, and absolute magnitude H. A more rehned approach would be 
to remove the known population from the population model. 

Mapping fpop to the best htting orbits in the (p, p) plane is quite com¬ 
plicated. Therefore, to compute such a prior we iteratively apply Bayesian 
inversion theory, i.e., the prior is computed as posterior of another prior 

f' ■ : 

J prior 


fprior {Pi p) = fpostiP, P) « f'prioriP^ P)fpop{q{p, P), e{p, p),i{p, p), H{p, /?)) . 


At this point we select fprior using simple geometric considerations. By as¬ 
suming that the Gartesian space of position and velocity is uniformly hlled, 
we have that p^ and p are uniformly distributed. Hence, fprioriP^ P) P^- 
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3.3. Uniform prior 


A common cho ice for f prior is a unifo rm di stribution. For iii s tance , such 


a prior is used by 


Virtanen et ah 


fl2nnih and 


Oszkiewicz et ah 


(120091). In¬ 


terestingly, it is possible to approximately obtain a prior distribution from 
the population model if we only consider the size distribution fsize{,H), i.e., 

fprior P fsizei,H'). 


For near-Earth asteroids the cumulative size distribution is expected to 


follow a power law 

(H) oc 

.0^^. with 77 somewhere between 0.35 and 

0.47 ( 

Bottke et ah. 

2002 ; 

Harris. 

2002 ; 

Stuart and Binzel. 

2004 

). For a 

given 

visual magnitude V, we have that fsizeiH) cx p { 

Bowell et ah. 

1989 

) and 


therefore fprior = According to the range of possible values for rj, this 

prior corresponds to a power law with exponent between —0.35 and 0.25, 
which is compatible with a uniform distribution. 


Though a uniform distribution might appear overly simplistic, we have 
found that this it does a good job of identifying potential impactors (e.g., see 
Sec. [6|) and distinguishing interesting objects from more common ones, such 
as main belt asteroids. 


4. Impact probability computation 

For each grid point, as a result of the constrained least squares fit to the 
observations, we obtain not only a best fit value Aij of the attributable, but 
also the corresponding uncertainty expressed by a covariance matrix 
This uncertainty is easily mapped to the orbital element space and therefore 
allows us to compute the probability Pij of an Earth impact for that orbit. 
Since for a fixed grid point the orbital uncertainty is small and the horizon 
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for impacts is of th e order of we e ks, we compute the impact pr obability with 


a linear approach fiMilani et ah 


2002 


Farnocchia et al 


r 


2015h 


Given fpost we can assign each grid point a weight Wij oc fpost{Pi, Pj), 
where — 1- Then, the impact probability can be easily computed as 


P = 


5. Generation of Monte Carlo orbital samples 

To better visualize the orbital distribution, it is convenient to generate 
Monte Carlo orbital samples distributed according to the probability distri¬ 
bution computed above. The following procedure can be used to achieve this 
goal: 

(a) Randomly draw according to a uniform distribution between the min¬ 

imum and maximum values of p on the grid; 

(b) Randomly draw pk according to a uniform distribution between the min¬ 

imum and maximum values of p on the grid; 

(c) Randomly draw a value 7 according to a uniform distribution between 

0 and the maximum of fpost on the grid; 

(d) Compute a best-£t attributable Ak for (p, p) = {pk,Pk)', 

(e) If fpostiPk, Pk) > 7 go to (f), otherwise start again from (a); 

(f) Add normally distributed noise to Ak according to the covariance r_ 4 j, 

and add {Ak,Pk,Pk) to the sequence of orbital samples. 


11 








Normalized astrometric RMS Normalized astrometric RMS 




Figure 2: Normalized RMS as a function of {p,p) for 2008 TC 3 . Left panel is for the first 
tracklet (four observations over 43 min), right panel is for the first two tracklets (total 
of seven observations over 99 min). Astrometric weights are at 0.5". The dashed curve 
corresponds to parabolic heliocentric orbits and the dash-dotted line to grazing impacts. 

6. Examples 

We want to test the algorithm presented in the previous sections on real 
cases. On one hand it is important that systematic ranging provides an early 
detection of impacting asteroids. On the other hand we need to avoid gen¬ 
erating too many warnings for objects that might be very far from Earth. 
Thus, we show examples of impactors as well as main belt asteroids. Ta¬ 
ble [ 2 ] summarizes the results for the examples discussed in more detail in the 
following subsections. 

6.1. 2008 TC3 

As already discussed in the introduction, 2008 TC 3 was the hrst aster¬ 
oid discovered prior to impact. This object was hrst spotted and promptly 
followed up by R. A. Kowalski at the Mt. Lemmon observatory. In the 
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20 hours between discovery and impact hundreds of astrometric observa¬ 
tions were reported to the Minor Planet Center. These observations allowed 
the calculation of a least squares orbit and the prediction of the impact. 
How would systematic ranging have performed with the hrst observations of 
2008 TC 3 ? 

The left panel of Fig. |2] shows the normalized RMS, i.e., the square root 
of the target function Q in Eq. ([1]), as a function of {p, p) by using the hrst 
tracklet of four observations obtained by Kowalski. The nonlinearity of the 
orbit determination problem is evident from the level curves of the RMS, 
which clearly suggest that 2008 TC 3 could be on an impact trajectory. As a 
matter of fact, the corresponding impact probability, computed as described 
in Sec. IU is already signihcant with the hrst four observations (Table |2]). 
With our defaults data weights at 0.5", the impact probability is few percents 
when using a uniform prior or one based on the population model, and ~60% 
when using Jehreys’ prior. With more conservative weights at 1" Jehreys’ 
prior leads to a ~70% impact probability, while the uniform and population 
prior result in a probability of few in a thousand. 

With the addition of Kowalski’s second tracklet, three additional observa¬ 
tions, the orbit becomes much better constrained by the observation residuals 
(see right panel of Fig. [2]). The impact is now almost certain, especially with 
0.5" weights. 

By using the hrst two tracklets of data and adopting a uniform prior, 
we also generated Monte Carlo samples to hnd potential impact locations, 
as shown in Fig. |3l As expectable, the region of possible impact loca¬ 
tion contains the Sudan desert where 2008 TC 3 actually ended its journey 
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Figure 3: Potential impact locations for 2008 TC 3 based on the first two tracklets of data. 


(iJenniskens et al.l . 1200911 . 


6.2. 2014 A A 

Though 2014 AA was not extensively followed up after discovery, there 
are several similarities with 2008 TC 3 . In particular, 2014 AA was also dis¬ 
covered about 20 hours before impact by R. A. Kowalski at the Mt. Lemmon 
observatory. 

Figure m shows the normalized RMS for 2014 AA as a function of (p, p) 
by using the first tracklet of data (three observations) or the full available 
dataset (seven observations). Again, the RMS curves clearly indicate that an 
impact is possible. Systematic ranging immediately recognizes the potential 
impact from the hrst tracklet and make the impact almost certain when the 
second tracklet is added. Table |2] reports the computed impact probabilities, 
which are similar to those computed for 2008 TC 3 ,. 

The impact locations (see Fig. [S]) found by using the seven available ob¬ 
servations and adopting a uniform prior are statistically consistent with the 
detection by the infrasound component of the International Monitoring Sys- 
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Normalized astrometric RMS Normalized astrometric RMS 




Figure 4: Normalized RMS as a function of {p,p) for 2014 AA. Left panel is for the first 
tracklet (three observations over 28 min), right panel is for the first two tracklets (total 
of seven observations over 69 min). Astrometric weights are at 0.5". The dashed curve 
corresponds to parabolic heliocentric orbits and the dash-dotted line to grazing impacts. 



Figure 5: Potential impact locations for 2014 AA based on all available data, i.e., seven 
observations. 
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Figure 6: Normalized RMS as a function of {p,p) for Hayabusa. Left panel is for the first 
tracklet (three observations over 7 min), right panel is for the first two tracklets (total 
of seven observations over 34 min). Astrometric weights are at 1.0". The dashed curve 
corresponds to parabolic heliocentric orbits and the dash-dotted line to grazing impacts. 


tern (IMS) operated by the Comprehensive Nuclear Test Ban Treaty Orga¬ 
nization (CTBTO), which places the impact location in the Atlantic Ocean 


at a l atitude of abo ut 15° N and a longitude of about 43° W fICheslev et al. 


2015 


Brown], 


2Qli). 


6.3. Hayabusa & Hayabusa 2 

Because of the low number of actual impactors we can use as a validation 


we te sted our algorithm on Hayabusa and Hayabusa 2 flKawaguchi et al. 


20081 ) . 


le reentry of Hayabusa took place on 2010 June 13. 


Yamaguchi et al. 


fl201ll ) report ground-based observations acquired during the reentry phase, 


which ended with a landing in the South Australia outback. We used these 
observations for applying systematic ranging technique to Hayabusa. Fig¬ 
ure El shows the normalized RMS as a function of (p, p) by using the hrst 
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tracklet of data (three observations) or the full available dataset (six obser¬ 
vations). Table [2] gives the impact probability as a function of the number of 
observations. With a uniform prior the initial impact probability is 2.7 x 10“^ 
and increases to 17% with the second tracklet. With the population-based 
prior the initial impact probability is 2.0 x 10“^ and increases to 2% with 
the second tracklet. On the other hand, with the use of Jefrreys’ prior we 
immediately have a 94% impact probability that decreases to 24% with the 
second tracklet. This counterintuitive behavior is due to the presence of two 
local RMS minima with the hrst tracklet, one of which clearly corresponds 
to an impact and is ruled out by the second tracklet. 


Hayabusa 2 was launched on 2014 December After launch Hayabusa 2 
was initially put on a low Earth parking orbit, until a maneuver moved it onto 
a ballistic trajectory leaving the Earth. On its way out, on 2014 December 6, 
Hayabusa 2 happend to be in the held of view of the Pan-STARRS PSl survey 
(IWainscoat et al.L 1201411 . which collected three observations of the spacecraft 
over 43 min. The corresponding astrometry was submitted to the Minor 
Planet Center, which, since the observations were not recognized as belonging 
to the Hayabusa 2 spacecraft, put the object on the NEO Conhrmation PageL 
with the temporary designation PlOgJfF. The following day, 2014 December 
7, Pan-STARRS obtained three additional observations from which it became 
clear that the object was Hayabusa 2. Table |2]shows how systematic ranging 
recognizes Hayabusa 2 as potential impactor from the hrst three observations 
with an impact probability probability of few to ten percents, depending on 


®http://global.jaxa.jp/press/2014/12/20141203_h2af26.html 

®http://www.ininorplanetcenter.net/iau/NEO/ToConfirin.htinl 
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the data weights and the prior used. Then, the subsequent observations 
better constrain the orbit and lower the impact probability since the ballistic 
trajectory does not correspond to an impact but to the post-injection escape 
orbit. 


64. 2004 ASi 


Asteroid 2004 ASi was discovered by LINEAR ( Stokes et al.l . [ 2 OOO) on 


2004 January 13 flBlvthe et al.l. l2004a|l . The four discovery observations span 
71 min and show an unusually high curvature thus suggesting that 2004 ASi 
was close to the Earth and that an impact was possible. Follow-up observa¬ 
tions later showed that the object was much farther than originally thought, 
thus ruling out any impact, and that the observed curvature was due to large 
astrometric errors in the discovery observations. 


Virtanen and MuinonenI (120061) employed the statistical ranging technique 


(jVirtanen et al.l. 1200111 to analyze the 2004 ASi case. They found an 18% 


impact probability with 1.0" weights and concluded that such a large im¬ 
pact probability was mostly due to the poor quality of the discovery night 
observations. 

The left panel of Fig. [7] shows the application of systematic ranging to 
2004 ASi. The astrometric observations are weighted at 1.0". The RMS 
favors an impact solution while the true solution (marked with a cross) cor¬ 
respond to larger astrometric errors. The impact probabilities computed for 
the different choices of the prior distribution are of the order of 10%, see 
Table |2l Decreasing the weights at 1.5" gives signihcantly lower probabilities 
for the uniform and population-based priors, while with Jefrreys’ prior we 
still have a 12% impact probability. 
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Figure 7: Normalized RMS as a function of (p, p) for 2004 ASi (left panel) and 2004 FU 162 
(right panel). Astrometric weights are at 1.0". The dashed curve corresponds to parabolic 
heliocentric orbits and the dash-dotted line to grazing impacts. For 2004 ASi the cross 
marks the true solution. 

As shown by Table [3l we computed the p-value of the true solution, i.e., 
the fraction of less likely orbital solutions, for the different choices of the 
prior and different data weights. Jeffreys’ prior produces very small p-values 
thus indicating that the true solution is deemed as extremely unlikely by the 
computed probability distribution. For the uniform distribution we have 3% 
p-value with 1" weights, which is due to the large errors in the astrometric 
observations. As a matter of fact, relaxing the weights increases the p-value 
to 22%. On the other hand, the population-based prior compensates for 
the large astrometric errors since impacting solutions with low RMS are less 
likely in the population model than non-impacting solutions with a higher 
RMS. The p-values are 13.6% and 38.7% depending on the data weights. 

The conclusion fpr 2004 ASi is that the discovery tracklet, weighted at 
1 ", yields impact probabilities ~ 10% even as the trueth is compatible with 
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the data at the 3-14% level. This outcome is a manifestation of the large 
astrometric errors in the data. 


6.5. 2004 FUi62 

Asteroid 2004 FU ifi-;i was also disco vered by LINEAR on 2004 March 31 


f Blvthe et ah 


2004 


CheslevI (120051 1 hrst analyzed this object and esti¬ 


mated that there was a 0.3% probability of an impact and that the object 
was likely to pass within 10 Earth radii. 

The right panel of Fig. [7]shows our implementation of systematic ranging, 
while Table [2] reports the impact probabilities for different weights and prior 
distributions. Again, Jeffreys’ prior results in very high impact probabilities 
while the population-based prior gives the lowest impact chances. For 1" 
weights the impact probability is 0.7% - 0.9% with the uniform and the 


population based priors. An Earth encounter within 1 
with all three priors, thus conhrming the Endings of 


0 Earth radii is assured 


CheslevI ( 200511 . 


6.6. 2015 CV & 2015 BUg2 

The previous examples show how systematic ranging is capable of quickly 
identifying actual and potential impactors. However, there is the risk of 
having significant impact probabilities ( 10 “^ or higher) for a large number of 
objects, which would make it difficult to distinguish interesting objects from 
less remarkable ones. Here we show the example of two main belt objects 
that were recently posted on the NEC Confirmation Page when the initial 
observations were compatible with an NEC orbit: 2015 CV and 2015 BU 92 . 

2015 CV was discovered by the Catalina Sky survey, Mt. Lemmon sta¬ 
tion, on 2015 February 10, while 2015 BU 92 was discovered by Spacewatch 
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Normalized astrometric RMS 



Topocentric range [au] 


Figure 8 : Normalized RMS as a function of (p, p) for 2015 CV (left panel) and 2015 BU 92 
(right panel). Astrometric weights are at 0.5". The dashed curve corresponds to parabolic 
heliocentric orbits and the dash-dotted line to grazing impacts. Crosses mark the true 
solutions. 


(iMcMillanl . 1200011 . Kitt Peak observatory, on 2015 January 16 (IBressi and Tubbiolol . 


20151). Figure IH] shows the application of systematic ranging to these two 


asteroids. For both objects the region with reasonable astrometric errors 
extends from impacting solutions to more than 1 au from the observer. 

As shown by Table I2| and Table 13], in both cases Jeffreys’ prior provides 
signihcant impact probabilities and extremely low p-values for the true so¬ 
lutions. In fact, Jeffreys’ prior is proportional to the partials of the astro- 
metric residuals and the residuals do not change much for large topocentric 
distances while they are very sensitive when the object is close to the ob¬ 
server. Therefore, Jeffreys’ prior tends to favor small topocentric distances 
even if the astrometric errors are larger. On the other hand the uniform and 
the population-based priors give much more reasonable probabilities and p- 
values. 
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7. Discussion and future work 


Small asteroids are expected to s trike the Earth on a regular basis. For 


instance, according to 


Brown et al. 


fj2002[) impacts with 10 m sized objects 


take place every few years. Since such small asteroids are likely to be ob¬ 
served, if at all, only shortly before impact, it is important to be capable of 
an early recognition of hazardous objects. 

The shortness of the observation arcs typically available for these ob¬ 
jects leads to critical degeneracies in the orbit estimation proc ess and limits 


the a pplicability of the standard impact monitoring systems flMilani et ah 


20051) . The systematic ranging technique presented in this paper relies on 


Bayesian inversion theory, which overcomes these degeneracies and results in 
a posterior distribution on the orbit phase space. 

A critical point is the choice of a prior distribution for the poorly con¬ 
strained space of topocentric range and range rate. This choice is far from 
trivial and indeed is somewhat arbitrary. We analyzed three options: Jef¬ 
freys’ prior, a uniform prior, and a prior distribution based on the asteroid 
population model. 

Jeffreys’ prior does a good job at assigning high impact probabilities to 
actual impactors. However, this prior produces high impact probabilities 
also for less remarkable objects such as main belt asteroids. The reason for 
this behavior is that Jeffreys’ prior favors regions in the phase space that are 
closer to the observer. As a matter of fact, Jeffreys’ prior is proportional to 
the partials of the residuals with respect to topocentric range and range rate, 
which are large for small topocentric distances and vanish at large distances. 
Moreover, in the non-impacting examples we analyzed we found that the true 
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solution was often deemed as extremely unlikely by the posterior distribution 
obtained starting from Jeffreys’ prior. Therefore, we avoid the use of Jeffreys’ 
prior for systematic ranging. 


The prior based on the asteroid orbit and size distribution represents 
the highest hdelity and most complex choice. However, it is worth pointing 
out that a population model is necessarily limited by the hnite number of 
samples and is therefore affected by Poisson statistics errors. In particular, 
small impactors might have peculiar orbits that are deemed unlikely by the 
model, which could be too much of a risk when the main goal is to identify 
potential impactors. 

A uniform distribution represents a compromise between simplicity and 
reliability. Though the corresponding posterior distribution and impact prob¬ 
abilities are not as rigorous as the one resulting from the population-based 
prior, we found that a uniform prior is effective at flagging potential im¬ 
pactors and discriminating between hazardous asteroids and more distant 
ones, which is the most important task to be achieved by systematic rang¬ 
ing. Another advantage of the uniform prior is that it has little prejudice on 
the actual orbit. 

To complement the standard impact monitoring performed by JPL’s Sen- 
tr}^ and the University of Pisa’s NEODysl^ we have implemented an auto¬ 
mated system to continually analyze the objects on the Minor Planet Center’s 
NEO Conhrmation Page with systematic ranging. This system quickly com¬ 
putes an impact probability to identify potential impactors and notihes us 


^http://neo.jpl.nasa.gov/risk 
®http: / / newton.dm.unipi.it / neodys 
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of interesting results. In such cases, we have communicated with discoverers 
and other observers to solicitate a quality control of the reported astrometry 
and for acquiring follow-up observations, which are necessary to decrease the 
orbital uncertainties and conclusively determine if the impact is going to take 
place. The examples of 2008 TC 3 and 2014 AA, i.e., the only two asteroids 
discovered before striking the Earth, suggest that objects with an impact 
probability above 10“^ should be prioritized for prompt follow-up. In recent 
experience this happens a few to several times a month. 

Among the possible additional applications of the systematic ranging 
technique presented in this paper, one could be the identihcation poten¬ 
tial radar targets or mission-accessib le targets, e.g., for the NASA’s Asteroid 


Redirect Mission flAbell et al 


201411 . 
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Object 

Number 

Arc 

Weights 

Impact probability 


of obs 

length 


Jeffreys 

Uniform 

Population 

2008 TCa 

4 

43 min 

0.5" 

6.3 X 10"^ 

4.4 X 10-2 

2.4 X 10-2 

2008 TCs 

4 

43 min 

1.0" 

6.7 X 10"^ 

3.7 X 10-3 

3.5 X 10-3 

2008 TCa 

7 

99 min 

0.5" 

1.0 

1.0 

1.0 

2008 TCs 

7 

99 min 

1.0" 

9.7 X 10-^ 

9.8 X 10-1 

9.6 X 10-1 

2014 AA 

3 

28 min 

0.5" 

8.1 X 10"^ 

2.9 X 10-2 

3.1 X 10-2 

2014 AA 

3 

28 min 

1.0" 

8.4 X 10"^ 

3.3 X 10-3 

8.4 X 10-3 

2014 AA 

7 

69 min 

0.5" 

1.0 

1.0 

1.0 

2014 AA 

7 

69 min 

1.0" 

9.9 X 10"^ 

9.6 X 10-1 

9.0 X 10-1 

Hayabusa 

3 

7 min 

1.0" 

9.4 X 10"^ 

2.7 X 10-3 

2.0 X 10-3 

Hayabusa 

6 

34 min 

1.0" 

2.4 X 10"^ 

1.7 X 10-1 

2.2 X 10-2 

Hayabusa 2 

3 

43 min 

0.2" 

2.8 X 10"^ 

3.6 X 10-1 

4.2 X 10-1 

Hayabusa 2 

3 

43 min 

0.5" 

3.6 X 10"^ 

3.3 X 10-2 

7.2 X 10-2 

Hayabusa 2 

6 

26 h 

0.2" 

2.4 X 10“^3 

2.3 X 10-13 

2.2 X 10-13 

Hayabusa 2 

6 

26 h 

0.5" 

1.9 X 10-3 

1.9 X 10-3 

1.8 X 10-3 

2004 ASi 

4 

71 min 

1.0" 

2.2 X 10-^ 

1.2 X 10-1 

7.8 X 10-2 

2004 ASi 

4 

71 min 

1.5" 

1.2 X 10-^ 

1.3 X 10-2 

1.1 X 10-2 

2004 FUi62 

4 

44 min 

1.0" 

3.6 X 10-^ 

7.3 X 10-3 

8.7 X 10-^ 

2004 FUig2 

4 

44 min 

1.5" 

5.0 X 10-^ 

1.3 X 10-2 

1.0 X 10-3 

2015 CV 

4 

34 min 

0.5" 

8.5 X 10-3 

4.7 X 10-® 

9.9 X lO-"^ 

2015 CV 

4 

34 min 

1.0" 

1.7 X 10-^ 

1.5 X 10-^ 

4.6 X 10-3 

2015 BU92 

3 

36 min 

0.5" 

1.7 X 10-^ 

7.2 X 10-3 

3.5 X 10-1° 

2015 BU92 

3 

36 min 

1.0" 

4.8 X 10-^ 

4.0 X 10-^ 

8.5 X 10-3 


Table 2: Impact probabilities for the examples considered in Sec. |6] for different choices of 
fprior, different number of tracklets, and different data weights. 
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Object 

Weights 

Jeffreys 

p-value 

Uniform 

Population 

2004 ASi 

1 . 0 " 

3.1 X 10"® 

2.6 X 10-2 

1.4 X 10-1 

2004 ASi 

1.5" 

4.0 X 10"^ 

2.2 X 10-^ 

3.9 X 10-1 

2015 CV 

0.5" 

5.0 X 10"® 

4.9 X 10-^ 

4.2 X 10-1 

2015 CV 

1 . 0 " 

1.3 X 10-5 

4.9 X 10-1 

4.2 X 10-1 

2015 BU 92 

0.5" 

1.1 X 10-5 

9.1 X 10-1 

8.8 X 10-1 

2015 BU 92 

1 . 0 " 

2.3 X 10-5 

9.1 X 10-1 

8.9 X 10-1 


Table 3: p-values for some of the examples considered in Sec. |6] for different choices of 
fprior and different data weights. 
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